Large polaron formation induced by Rashba spin-orbit coupling 
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Here the electron-phonon Holstein model with Rashba spin-orbit interaction is studied for a two 
dimensional square lattice in the adiabatic limit. It is demonstrated that a delocalized electron at 
zero spin-orbit coupling localizes into a large polaron state as soon as the Rashba term is nonzero. 
This spin-orbit induced polaron state has localization length inversely proportional to the Rashba 
coupling 7, and it dominates a wide region of the 7- A phase diagram, where A is the electron-phonon 
interaction. 
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I. INTRODUCTION 



Spin manipulation and control is at the core of spin- 
tronics, a technology that uses the spin of the elec- 
trons, rather than their charge, to transfer and/or pro- 
cess information.^ The Rashba spin-orbit (SO) cou- 
pling arising in materials lacking structural inversion 
symmetry^ plays a leading role in this field because its 
strength can be tuned by an applied electric field and by 
specific material engineering methods. The SO induced 
lifted spin degeneracy may then be used in spin filtering 
devices and spin transistors. 

Whether the main effect of SO coupling is limited to 
the spin splitting or it is accompanied by substantial 
modifications in other electronic properties, which could 
be detrimental for the spin propagation, is of course cru- 
cial for the functioning of spin-based devices. In this 
respect, an important issue calls into play the role of the 
SO interaction on the coupling of electrons to the lattice 
vibrations (phonons). In particular, a sensible problem is 
whether the polaron, that is the quasiparticlc composed 
by the electron and its phonon cloud, is strengthened or 
weakened by the Rashba SO interaction. 

In previous works, an enhancement of the polaronic 
character has been obtained for a two-dimensional (2D) 
electron gas with linear Rashba coupling for both short- 
range (Holstein model Ref.Q) and long-range (Frohlich 
model Ref.[f|) electron-phonon (el-ph) interactions^— 
On the contrary, a recent calculation on the 2D tight- 
binding Holstein-Rashba model on the square lattice has 
shown that a large el-ph interaction gets effectively sup- 
pressed by the Rashba SO coupling £ At present there- 
fore the role of the Rashba SO coupling on the polaron 
properties is not clear, and different models and approx- 
imations appear to give quite contradicting results. 

In this article the tight-binding Holstein-Rashba model 
for one electron coupled to adiabatic phonons is consid- 
ered and the corresponding non-linear Schrodingcr equa- 
tion for the polaron wave function is solved numerically. 
It is shown that, for el-ph couplings such that the elec- 
tron is delocalized in the zero SO limit, the Rashba term 
creates a large polaron state, with polaron localization 
length inversely proportional to the SO strength. Fur- 
thermore, the small polaron regime appearing at large 



el-ph couplings and zero SO gets weakened (or even sup- 
pressed) for sufficiently strong SO couplings. Hence, the 
Holstein-Rashba polaron is strengthened or weakened by 
the SO interaction depending on whether the el-ph cou- 
pling is respectively weak or strong, thereby reconciling 
the different trends reported in Rcfs.[6.9] into one single 
picture. 



II. MODEL 



By presenting the spinor operator ^ R 



where ct>„, creates an electron with spin a 
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R, the tight-binding Holstein-Rashba Hamiltonian on the 
square lattice can be written as H — Hq + H p h + H e i^ p h, 
where 



Hn = - 



R 

*2 H (*R°V*H+* - *k^*R+y 



H.c. 



(1) 



is the lattice Hamiltonian for a free electron with trans- 
fer integral t and SO coupling 7. a x and a y are Pauli 
matrices. The lattice constant is taken to be unity, 
and x and y are unit vectors along the x and y direc- 
tions, respectively. The Hamiltonian ([1]) is easily di- 
agonalized in momentum space, and the resulting elec- 
tron dispersion is composed of two branches: E^ = 
— 2t[cos(k x ) + cos(k y )] ±7y / sin(fc r ) 2 + sin(£; y ) 2 . The low- 
est branch, Er, has a four- fold degenerate minimum 
E = — 4^1 + 7 2 /(8t 2 ) for momenta k = (±fc ,±fc ) 

with ko = arctan[7/(\/8i)]-— The Hamiltonian for Ein- 
stein phonons with mass M and frequency loq is given 
by: 
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where Pr and are impulse and displacement phonon 
operators. Finally, the el-ph Hamiltonian contribution is 
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FIG. 1: (Color online) Total energy difference AE — E — Eq for the adiabatic Holstein-Rashba model as a function of the el-ph 
coupling A and for different values of the SO Rashba interaction 7. Eo is the ground state energy for A = 0. Different symbols 
refer to different solutions of the nonlinear Schrodinger equation ([5]), and the ground state is given by the solution with lower 
AE values. Insets: corresponding electron density probability at R = 0. 



where g is the el-ph interaction matrix element. 

The (quasi-) 2D materials and heterostructures which 
display non-zero Rashba couplings (semiconductor quan- 
tum wells, surface states of metals and semimetals) are 
wide electron bandwidth systems with t of the order of 1 
eV, while the typical phonon energy scale is of the order 
of few to tens meVJ^- These systems are expected there- 
fore to be well within the adiabatic regime u>o/t <C 1. 
In the following, however, only the strict adiabatic limit 
(jJo/t — is considered, which simplifies considerably the 
problem and, as shown below, permits to identify the crit- 
ical parameters governing the electron localization tran- 
sitions. 

The adiabatic limit uio/t = is obtained formally 
from Eqs. ([2]) and ([3]) by setting M — > 00 and keep- 
ing K = Muiq finite. Since for M — > 00 the phonon 
kinetic energy is zero, the ground state in the adiabatic 
limit is obtained by finding the displacement configura- 
tion which minimizes the total energy E = (H), 
where the brackets mean the expectation value with re- 
spect to the electron wave function and the lattice dis- 
placement. Hence, since by Hellmann-Feynman theorem 

X^ = V2AM) ff( l 0|*R*R|V') /K, the ground state en- 
ergy becomes 

E GS = (1>\H \1>) - E P ^<V>|*k* R |V> 2 , (4) 

R 

where Ep = g 2 /uiQ is independent of M, and \ip) = 
Sr q ^Ra c Ra|0)- The ground state electron wave func- 
tion 4>ji a can be found from Eq.Q by applying the 
variational principle, leading to the following non-linear 
Schrodinger equation: 

£$R = - t ( $ R+«* + *R+«y) - 2E P I $ R I 2 $ R 

n=± 

n(a y $ R+n z - a x $ R+n y), (5) 

1 n=± 



where $ R = ((f>* RV R; )+ and e = E GS + E P £ R |$ R | 4 . 
Finally, the ground state energy Eqs is obtained by solv- 
ing Eq.© iteratively, with J2 R a |0r q | 2 = 1, and by in- 
serting the resulting wave function into Eq.@. 



III. RESULTS 

Solutions of ([5]) for lattices of N = 101 x 101 sites 
are plotted in Fig. [1] BjS ci function of the el-ph cou- 
pling constant A = Ep/{At) = g 2 /(Atwo) and for four 
different values of 7. For 7 = 0, Fig. [TJa), we re- 
cover the well-known behavior of the adiabatic Holstein 
model in two-dimcnsions?i2, a dclocalized solution with 
Eqs — Eq = —At (filled circles) extending to the whole 
range of A values considered, and a localized one (filled 
squares) having energy lower than Eg for A > A* = 0.835. 
The delocalized/localized nature of the solutions is illus- 
trated in the inset of Fig.[TJa) where the electron density 
probability |<&r| 2 = J2 a I^Ra| 2 is plotted for R = 0. 
The solution having lower energy for A > A* is a small 
polaron state, with more than 90% of its wave function 
localized at the origin. 

Let us now consider the 7 > case. As shown in 
Figs.QJb)-(c), a nonzero Rashba term gives rise to a new 
feature absent for 7 = 0. Namely, besides the two solu- 
tions already discussed for the 7 = third solution 
appears (filled triangles), which has lower energy than 
the delocalized and small polaron states in a region of 
intermediate values of A. It is thus possible to identify a 
second critical coupling, A**, such that for A** < A < A* 
the ground state is given by this third solution. Further- 
more, the transition to the small polaron state (identi- 
fied by A*) gets shifted to larger el-ph couplings as j/t 
increases, thereby confirming the results of Ref.Q ob- 
tained by a different method and for Uo/t ^ 0. A map of 
the behavior of A* and A** as 7 is varied is reported in 
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FIG. 2: (Color online) Phase diagram of the 2D adiabatic 
Holstein-Rashba model. The A** and A* transition lines are 
the phase boundaries separating the different states of the 
polaron. The dashed curve has been obtained from the max- 
imum of d 2 E/X 2 and identifies a smooth crossover from large 
to small polaron for large y/t values. The solid line is the 
variational result of Eq. (|12l) . The graded gray (violet) scale 
refers to the polaron density probability at R = 0. 



the y/t- A phase diagram of Fig. [21 where the filled circles 
are the calculated values of A**, while the filled squares 
mark the onset of the small polaron regime (A*)J^ The 
resulting diagram is therefore composed of three sepa- 
rate regions: a delocalized electron with Eqs = Eq for 
y/t > A** (white region), a small polaron state for large 
el-ph couplings (A > A*), and a new ground state in the 
region comprised between the A** and A* lines. 

As it can be inferred from the insets of Fig. [T] and from 
the gray (violet) scale of Fig. [21 in this region the den- 
sity probability at R = 0, |$o| 2 ; is lower than the small 
polaron solution, but substantially larger than zero as 
long as 7 7^ 0, and increasing with y/t. The region be- 
tween the A** and A* lines identifies therefore a large po- 
laron state created by the SO interaction, with a localized 
wave function which may extent over several lattice sites. 
The large polaron nature of this solution is substantiated 
in Fig. [31 where the polaron localization radius Rp, ex- 
tracted from a fit of |$r| 2 to exp(— |R|/i?p) (see inset), 
is plotted as a function of y/t for A = 0.4, 0.6, and 0.8. 
Although a numerical evaluation of Rp for y/t — > is 
hampered by the finite size of the lattice, Rp turns out to 
be approximately proportional to t/y, suggesting there- 
fore that the large polaron evolves continuously towards 
a delocalized electron as y/t — > 0. 

Further insight on the large polaron state, and in par- 
ticular on its behavior as y/t — > 0, can be gained by a 
simple variational calculation in the continuum. In fact, 
as long as Rp is much larger than the lattice constant 
(Rp 3> 1) then an upper bound for Eqs can be obtained 
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FIG. 3: (Color online) Polaron radius Rp of the large polaron 
state as a function of y/t and for different el-ph couplings A. 
Inset: Density probability (symbols) of the large polaron for 
A = 0.4 and y/t = 0.25, 0.5, and 1.0 as a function of distance 
R — |R| along the (1,0) direction. The solid lines are fits to 
exp(-|R|/i? P ). 



from a minimization of the energy functional 

£[$(r)] = Jdr&(r)[tp 2 +y{a y p x -<T x p y )]<P{r) 

-E P JdvMv)\\ (6) 

where p q = —id/dq is the electron momentum operator 
(q = x,y and H = 1) and p 2 = p 2 + p 2 . In the above 
expression, ^(r) is a suitable ansatz for the ground state 
spinor, which is assumed to vary slowly over distances 
comparable to the lattice spacing. In writing Eq.©, only 
the lowest order terms in the lattice constant have been 
retained, which amounts to consider a parabolic band 
with a Rashba coupling linear in the momentum oper- 
ators. One can then use for $(r) an ansatz which has 
been already introduced in studying the effects of a lin- 
ear Rashba term on the 2D Frohlich polaron and the 2D 
hydrogen atom&il 



$(r) = Acxp(— ar) 



Mbr) 
Ji(6r) e i 



(7) 



Here, r = |r| and ip is the azimuthal angle, A is a nor- 
malization constant, Jo an d J\ are Bessel functions, and 
a and b are variational parameters. By using ([7|) and the 
properties of the Bessel functions, equation ^ reduces 
to 



£ = t(a 2 + b 2 )-yb- P Jo 



-4ar 



F(br) 



t(a 2 



yb 



2- [J- dr 

2E P a 2 , 
m 



-2ar 



F(br)\ 



(8) 



where F(br) = Jo(br) 2 + J\(br) 2 . The second equality 
stems from assuming a <C b, which is the relevant limit of 
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the large polaron regime. Minimization of £ with respect 
to a and b leads to two possible solutions: 6 = 7/ (2t) and 
a = 0, which corresponds to a dclocalized electron with 
£min = E Q = -At - j 2 /(At), and 

( k \ , 7/(2*) 

fl = 6e xp (-1 - -j , 6 = — 4Aexp h2 _ ^ , 

(9) 

which represents the large polaron solution with 

ftnin - #o = ~y exp (-2 - ^) , (10) 

for A small. Since Eq. (|10[) is an upper bound for AE'qs = 
Eqs — Eq, then that the large polaron state has energy 
always lower than the delocalized electron. Furthermore, 
by realizing that the variational parameter a represents 
the polaron radius through a = l/(2Rp), it turns out 
from Eq.(|9j that Rp scales as t/j, in agreement therefore 
with the results of Fig. [3] 

The finding that a large polaron is formed for 7/t 7^ 
is in accord with the observation of Ref.0 that per- 
turbation theory breaks down in the adiabatic limit 
for any finite A. This breakdown basically stems from 
the one-dimensional-like divergence of the density of 
states (DOS) of a parabolic band with linear Rashba 
coupling 

Although the variational result presented above cor- 
rectly predicts the appearance of the large polaron state 
as soon as j/t 7^ 0, it fails nevertheless in describing the 
A** transition line of Fig. [2] separating the large polaron 
state from the delocalized solution. This is because the 
lowest order expansion in the lattice constant of Eq. (|6]) 
neglects higher order powers of the momentum operator 
arising from the lattice Rashba term, which shift the van 
Hove divergence of the DOS from Eq to higher energies^ 
thereby making the perturbation theory non-singular. To 
investigate this point within the variational method, it 
suffices to expand the discrete Hamiltonian up to the 
third order in the lattice constant. This corresponds to 
add to the energy functional (|6]) the following contribu- 
tion 

£'[*(R)] = I J dv&{v)(a x f y - a y p 3 x )$(r), (H) 

which, by using again the ansatz (JT)) and for a <C b, leads 
to the third order correction term £' = (7/8)(fe 3 + 3a 2 b — 
na 3 ) to Eq.©. It is then easy to shown that [£ +£'} min — 
E is negative [with E = -At - 7 2 /(4t) + 7 4 /(128t 3 )] as 
long as 7/t < A** r , where for A small 

A-=8/f exp(-l-^). (12) 

Although Eq. (fT2"]) provides only a lower bound for A** 
(solid line in Fig. [2]), it shows nevertheless that, as 7/t is 
enhanced for fixed A, the transition from the large po- 
laron to the dclocalized electron state originates from 
higher order of the SO interaction than the linear Rashba 
coupling. 



IV. DISCUSSION AND CONCLUSIONS 



Let us discuss now the significance of the results re- 
ported above for materials of interest and possible con- 
sequences for spintronics applications. First of all, it is 
important to identify the region in the phase diagram of 
Fig. [5] where realistic values of 7/t and A are expected 
to fall. This is easily done by realizing that the largest 
Rashba SO coupling to date is that found in the surface 
stats of Bi/Ag(lll) surface alloys^ for which j/t « 1.4 
can be estimated. Other 2D systems and heterostruc- 
tures have lower or much lower 7/t values. Concerning 
the coupling to the phonons, a survey^ on the el-ph in- 
teraction at metal surfaces evidences that A is usually 
lower than 0.6-0.7 (see also Ref.[I3|), at least for the sur- 
face states with large SO splittings (i.e. Ag, Cu, Bi). It 
is therefore a rather conserving assumption to confine to 
7/t < 1 and A < 1 the region of interest for the micro- 
scopic parameters which, as shown in Fig. [5J is substan- 
tially dominated by the SO induced large polaron state. 
Hence, upon tuning of the Rashba SO coupling, a delocal- 
ized electron at 7/t = can in principle be changed into a 
self-trapped large polaron state for 7/t > 0, with obvious 
consequences on the spin propagation in the system. In 
passing, it is worth noticing that the small polaron regime 
instead is affected rather weakly by the SO interaction 
for 7/t < 1, while its weakening gets pronounced only for 
unrealistically large values of 7/t (see also Fig.[T]). 

Before concluding, it is important to discuss a last im- 
portant point. Although the adiabatic limit employed 
here allows for a clear identification of the A* and A** 
transition lines, the energy gain associated to the large 
polaron formation becomes very small in the weak cou- 
pling and small SO limits [see Eq. (flU)) ]. In this regime, 
the inclusion of quantum fluctuations which arise as soon 
as ujo/t 7^ may wash out completely any signature (like 
e.g. an anomalous enhancement of the electron effective 
mass m*) of the large polaron state, even for u>o/t small, 
while they should remain visible for larger A and 7/t val- 
ues. For a more complete description of the SO effects on 
the Holstcin-Rashba polaron, it is therefore necessary to 
extend the study to the non-adiabatic regime uj /t 7^ 0, 
by keeping however in mind that, as discussed above, 
relevant materials have u)o/t <C 1. 

In summary, the complete phase diagram of the 2D 
adiabatic Holstein el-ph Hamiltonian in the presence of 
Rashba SO coupling has been calculated. It has been 
shown that a self-trapped large polaron state is created 
by the SO interaction in a wide region of the phase dia- 
gram, and that its localization radius can be modulated 
by the SO coupling. This result implies that, for realistic 
values of the microscopic parameters, the appearance of 
a self-trapped large polaron state is a potentially detri- 
mental factor for spin transport. 
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